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The causal properties of curved spacetime, which underpin our sense of time in 
gravitational theories, are defined by the null cones of the spacetime metric. In 
classical general relativity, it is assumed that these coincide with the light cones de- 
termined by the physical propagation of light rays. However, the quantum vacuum 
acts as a dispersive medium for the propagation of light, since vacuum polarisa- 
tion in QED induces interactions which effectively violate the strong equivalence 
principle (SEP). For low frequencies the phenomenon of gravitational birefringence 
occurs and indeed, for some metrics and polarisations, photons may acquire super- 
luminal phase velocities. In this article, we review some of the remarkable features 
of SEP violating superluminal propagation in curved spacetime and discuss recent 
progress on the issue of dispersion, explaining why it is the high-frequency limit of 
the phase velocity that determines the characteristics of the effective wave equation 
and thus the physical causal structure. 

1 Introduction 

In general relativity, the nature of time and causality are determined by the prop- 
erties of the null cones of the curved spacetime manifold. A physical realisation 
of these geometric null cones is provided in classical electrodynamics by the light 
cones traced by the propagation of light rays. As the title of this conference implies, 
'Time and Matter' are therefore intimately related. 

In quantum theory, however, the relationship is more subtle. Quantum effects 
such as vacuum polarisation in QED induce interactions which effectively violate 
the strong equivalence principle and cause the quantum vacuum in the presence of 
gravity to act as a dispersive medium for the propagation of light. At least for low 
frequencies, the phenomenon of gravitational birefringence occurs and moreover, for 
some metrics and polarisations, photons may acquire superluminal phase velocities. 
This forces a reassessment of the identification of light cones with null cones and 
raises the question of how causality and quantum theory can be reconciled in general 
relativity. 

Research into this effect began with a key paper of Drummond and Hathrell^ in 
which they calculated the one-loop vacuum polarisation contribution to the QED 
effective action in a background gravitational field. They found the following mod- 
ification to the free Maxwell action: 



where a,b,c,d are constants of 0{a) and m is the electron mass. The important fea- 
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ture is the direct coupling of the electromagnetic field to the curvature. This is an 
effective violation of the strong equivalence principle (SEP), which is the dynamical 
ansatz that the 'laws of physics' should be the same in the local incrtial frames at 
each point in spacetime. (The weak equivalence principle, viz. the existence of such 
local inertial frames, is in contrast a fundamental assumption underlying the struc- 
ture of general relativity.) More precisely, the SEP requires that elcctromagnctism 
is minimally coupled to gravity, i.e. through the connections only, independent of 
the curvature. The effective action Eq.(l) shows that while this principle may be 
consistently imposed at the classical level, it is necessarily violated in quantum 
electrodynamics. It is the quantum violation of the SEP that allows the physical 
light cones to differ from the geometric null cones. 

The new SEP- violating interactions in Eq.(l) affect the propagation of light 
and modify the physical light cones. Using the techniques of geometric optics, 
Drummond and Hathrell showed that the new light cones are given by 

e - ^R^^k^k^ + ^^ k^^k^a'^aP = (2) 

where k^ is the wave vector and the polarisation. Using the Einstein equa- 
tions, this can be expressed as the sum of a 'matter' and a purely 'gravitational' 
contribution as follows: 

- -^(26 + 4c)T^Afc''fc^ + ^C^^Xph^k^a^aP = (3) 

The first term involves the projection of the energy-momentum tensor which ap- 
pears in the weak energy condition; this contribution is universal, appearing in 
the modified light cone condition for light propagating subluminally in a variety 
of backgrounds such as classical electromagnetic fields or finite temperature 2.3.4,5_ 
The second term depends on the Weyl tensor and is uniquely gravitational; since it 
depends explicitly on the polarisation, the modified light cones exhibit gravitational 
birefringence. 

It is also instructive to express Eqs.(2),(3) in the Newman-Penrose formalism. 
Introducing a null tetrad with basis vectors n^, and together with 
the corresponding components of the Ricci and Weyl tensors <I>oo = —^R^v^-^f-", 
= — CuiyXpiP-m'^i^mP etc. (for details, see ref.'^), and choosing to coincide 
with the direction of propagation, i.e. A;^ = Lolfj,, we find 
,2 (46 + 8c)w2 4cw2 

This representation makes it clear that the contribution from the Weyl tensor 
changes sign for the two physical transverse polarisations a^. It follows imme- 
diately that for Ricci- flat spacetimes, both timelike and spacelike values of k'^ are 
possible. In other words, Eqs.(2),(3);(4) necessarily imply the existence of superlu- 
minal propagation. Physical photons no longer follow the geometrical null cones, 
but instead propagate on the effective, polarisation-dependent, light cones defined 
above. 

Many examples of the Drummond-Hathrell effect in a variety of gravitational 
wave, black hole and cosmological spacetimes have been studied i'^'^''^. The (Ricci- 
flat) black hole cases are particularly interesting, and it is found that for photons 
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propagating orbitally, the light cones are modified such that superluminal propa- 
gation occurs. For radial geodesies in Schwarzschild spacetime, however, and the 
corresponding principal null geodesies in Rcissncr-Nordstrom and Kerr, the light 
cones are unchanged. The reason is simply that if we choose the standard Newman- 
Penrose tetrad in which £^ is tangent to the principal null geodesic, the only non- 
vanishing component of the Wcyl tensor is ^'2 since these black hole spacetimcs arc 
all Petrov type D, whereas the modification to the light cone condition involves only 
^Q. Superluminal propagation is also predicted in the (Weyl-flat) FRW cosmologi- 
cal spacetimcs, with the correction to the speed of light increasing as towards 
the initial singularity. Another interesting case involves the Bondi-Sachs metric 
describing gravitational radiation from an isolated source, where the magnitude of 
the superluminal effect is related to the peeling theorem for the Wcyl tensor. 

The existence of superluminal propagation in QED in curved spacetime of course 
raises immediate questions as to the realisation of causality. The purpose of this 
article is to review our research programme on photon propagation in gravitational 
fields with particular emphasis on the issue of causality and the consistency of 
quantum field theory with classical gravitation. We begin, in section 2, by con- 
sidering carefully the implications of the Drummond-Hathrcll effective action as it 
stands, reviewing the bimetric interpretation of the light cones, the realisation of 
stable causality in general relativity, possible time machine constructions and the 
consequences for event horizons. 

However, the action Eq.(l) is only the lowest-order term in an expansion of the 
full effective action in powers of derivatives. That is, results derived from it are valid 
only in a low- frequency approximation A ^ Ac, where Ac is the electron Compton 
wavelength. The inclusion of terms of higher orders in derivatives shows that the 
propagation is in fact dispersive^. In section 3, we discuss the precise definition of 
the 'speed of light' and present a proof that the wavcfront velocity, which is the 
relevant speed of light for causality, can be identified as the high-frequency limit 
of the phase velocity, Vph{oo). This means that a resolution of the issues raised in 
section 2 concerning causality depends on an explicit calculation of the light cones 
for high-frequency propagation. In section 4, we present a recently derived extension 
of the effective action valid to third order in curvatures and field strengths and to 
all orders in derivatives^. The resulting light cone is derived and some potential 
special features are described. Finally, however, we argue on the basis of a detailed 
comparison with photon propagation in background magnetic fields that a further, 
non-perturbative contribution to the effective action may ultimately control the 
high-frequency limit and we close with an assessment of the prospects for a final 
resolution of the question of dispersion and causality for QED in curved spacetime. 



2 Causality and Superluminal Propagation 

In this section, we discuss the implications for causality of photon propagation based 
on the Drummond-Hathrell action Eq.(l) and its associated light cones, setting 
aside for the moment the issue of dispersion. 
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Figure 1. A superluminal {v > 1) signal OA which is forwards in time in frame S is backwards 
in time in a frame S' moving relative to S with speed u > ^. However, the return path with the 
same speed in <S arrives at B in the future light cone of O, independent of the frame. 

2.1 Superluminal Propagation in Special and General Relativity 

It is generally understood that superluminal propagation in special relativity leads 
to unacceptable violations of causality. Indeed the absence of tachyons is tradi- 
tionally employed as a constraint on fundamental theories. We therefore begin by 
reviewing some basic features of superluminal propagation in order to sharpen these 
ideas in preparation for our subsequent discussion of causality in general relativity. 

The first important observation is that given a superluminal signal we can always 
find a reference frame in which it is travelling backwards in time. This is illustrated 
in Fig. 1. Suppose we send a signal from O to A at speed u > 1 (in c = 1 units) 
in frame S with coordinates (t,x). In a frame S' moving with respect to S with 
velocity u > ^, the signal travels backwards in t' time, as follows immediately from 
the Lorentz transformation. ° 

The important point for our considerations is that this by itself does not nec- 
essarily imply a violation of causality. For this, we require that the signal can be 
returned from A to a point in the past light cone of O. However, if we return the 
signal from A to B with the same speed in frame S, then of course it arrives at 
B in the future cone of O. The situation is physically equivalent in the Lorentz 
boosted frame S' - the return signal travels forward in t' time and arrives at B 
in the future cone of O. This, unlike the assignment of spacetime coordinates, is a 
frame- independent statement. 

The problem with causality arises from the scenario illustrated in Fig. 2. Clearly, 

"^From the Lorentz transformations, we have t'^ = 7(n)tj4(l — uv) and x'^ = 'y{u)xA{^ ~ ^) ^'^^ 
the situation realised in Fig. 1, we require both x'^ > and t'^ < 0, that is ^ < u < v, which 
admits a solution only if > 1. 
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Figure 2. A superluminal (v > 1) signal OC which is backwards in time in frame S is returned at 
the same speed to point D in the past light cone of O, creating a closed time loop. 

if a backwards-in-time signal OC is possible in frame S, then a return signal sent 
with the same speed will arrive at D in the past light cone of O creating a closed 
time loop OCDO. The crucial point is that local Lorentz invariance of the laws of 
motion implies that if a superluminal signal such as OA is possible, then so is one 
of type OC, since it is just given by an appropriate Lorentz boost (as in Fig. 1). 
The existence of global inertial frames then guarantees the existence of the return 
signal CD (in contrast to the situation in Fig. 1 viewed in the S' frame). 

The moral is that both conditions must be met in order to guarantee the occur- 
rence of unacceptable closed time loops - the existence of a superluminal signal and 
global Lorentz invariance. Of course, since global Lorentz invariance (the existence 
of global inertial frames) is the essential part of the structure of special relativity, 
we recover the conventional wisdom that in this theory, superluminal propagation 
is indeed in conflict with causality.'' 

The reason for presenting this elementary discussion is to emphasise that the 
situation is crucially different in general relativity. The weak equivalence principle, 
which we understand as the statement that a local inertial frame exists at each 
point in spacetime, implies that general relativity is formulated on a Riemannian 
manifold. However, local Lorentz invariance alone is not sufficient to establish the 
link between superluminal propagation and causality violation. This is usually 
established by adding a second, dynamical, assumption. The strong equivalence 
principle (SEP) states that the laws of physics should be identical in the local 
frames at different points in spacetime, and that they should reduce to their special 
relativistic forms at the origin of each local frame. It is the SEP which takes over 
the role of the existence of global inertial frames in special relativity in establishing 

''Notice that this is not in contradiction with the occurrence of superluminal propagation in flat 
spacetime with Casimir plates^"'^^'^^ , since the modification of the spacetime geometry by the 
plates removes global Lorentz invariance. 



Shore: submitted to World Scientific on February 7, 2008 5 



the incompatibility of superluminal propagation and causality. 

However, unlike the weak equivalence principle, which underpins the essential 
structure of general relativity, the SEP is merely a simplifying assumption about the 
dynamics of matter coupled to gravitational fields. Mathematically, it requires that 
matter or electromagnetism is minimally coupled to gravity, i.e. with interactions 
depending only on the connections but not the local curvature. This ensures that 
at the origin of a local frame, where the connections may be Lorentz transformed 
locally to zero, the dynamical equations recover their special relativistic form. In 
particular, the SEP is violated by interactions which explicitly involve the curvature, 
such as those occurring in the Drummond-Hathrell action (1) and the consequent 
modified light cones. 

The question of whether this specific realisation of superluminal propagation 
is in conflict with causality is discussed in section 2.4 using the concept of stabele 
causality described in ref.^^. Notice though that by violating the SEP, we have 
evaded the necessary association of superluminal motion with causality violation 
that held in special relativity. Referring back to the figures, what is established 
is the existence of a signal of type OA, which as we saw, does not by itself imply 
problems with causality even though frames such as S' exist locally with respect 
to which motion is backwards in time. However, since the SEP is broken, even if 
a local frame exists in which the signal looks like OC, it does not follow that a 
return path CD is allowed. The signal propagation is fixed, determined locally by 
the spacetime curvature. 



2.2 Geometric Optics 

The most direct way to deduce the form of the light cones for QED in curved 
spacetime is to use geometric optics. This starts from the ansatz 

(^^ + ieBi, + ...) exp(i-) (5) 

in which the electromagnetic field is written as a slowly-varying amplitude and a 
rapidly-varying phase. The parameter e is introduced as a device to keep track of 
the relative order of magnitude of terms, and the Bianchi and Maxwell equations 
are solved order-by-order in e. 

The wave vector is identified as the gradient of the phase, fc^ = We also 
write = Aa^, where A represents the amplitude itself while specifies the 
polarisation, which satisfies k^a^ = 0. 

Solving the usual Maxwell equation D^F^'^ = 0, we find at 0(l/e), 

= (6) 

while at 0(1), 

k'^D^a'' = (7) 

and 

k''D^{\nA) = ~D^k'' (8) 
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Eq.(6) shows immediately that fc^ is a null vector. From its definition as a gradient, 
we also see 

k'^Df.k'' = k^D'^k^ = ^D'^e = (9) 

Light rays, or equivalently photon trajectories, are the integral curves of k^, i.e. the 
curves a;'*(s) where dx^/ds = k^^. These curves therefore satisfy 

This is the geodesic equation. We conclude that for the usual Maxwell theory 
in general relativity, light follows null geodesies. Eqs.(9),(7) show that both the 
wave vector and the polarisation are parallel transported along these null geodesic 
rays, while Eq.(8), whose r.h.s. is just (minus) the optical scalar 9, shows how the 
amplitude changes as the beam of rays focuses or diverges. 

2.3 Bimetricity 

The same method is applied to the modified Maxwell equation derived from the 
effective action Eq.(l): 



(11) 



from which the new light cone condition 



fc2 - ^R^^kf'k^ + ^R^^Xflk^'k^a'^aP = (12) 



follows immediately. Since this new light cone relation is still homogeneous and 
quadratic in k^, we can write it as 

Q^'^ki^k. = (13) 

defining Q^" as the appropriate function of the curvature and polarisation. 

Now notice that in the discussion of the free Maxwell theory, we did not need to 
distinguish between the photon momentum p^, i.e. the tangent vector to the light 
rays, and the wave vector k^ since they were simply related by raising the index 
using the spacctimc metric, = In the modified theory, however, there is 

an important distinction. The wave vector, defined as the gradient of the phase, is 
a covariant vector or 1-form, whereas the photon momentum/tangent vector to the 
rays is a true contravariant vector. The relation is non-trivial. In fact, given k^, 
we should define the corresponding 'momentum' as 



p'' = G'^'K (14) 

ds 



and the light rays as curves x'^{s) where = p^. This definition of momentum 



satisfies 

G^^p^p" = g^^k^k^ = (15) 

where G = defines a new effective metric which determines the light cones 
mapped out by the geometric optics light rays. (Indices are always raised or lowered 
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using the true metric g^v) The ray velocity v^ay corresponding to the momentum 
p^, which is the velocity with which the equal- phase surfaces advance, is given by 
(defining components in an orthonormal frame) 



\P\ _ d]g 
dt 



"""^y - 7^ ~ ~ ^^^^ 



along the ray. This is in general different from the phase velocity 

V,. = ^ (17) 

This shows that photon propagation for QED in curved spacetime can be char- 
acterised as a bimetric theory^ - the physical light cones are determined by the 
effective metric G^i^ while the geometric null cones are fixed by the spacetime met- 
ric g^u- 



2.4 Stable Causality 

The bimetric formulation is the most natural language in which to discuss whether 
the superluminal velocities predicted by the Drummond-Hathrell action are com- 
patible with our usual idea of causality. We have already seen that with SEP 
violation in general relativity, the arguments that in special relativity led to the 
incompatibility of superluminal propagation and causality are no longer valid. Su- 
perluminal motion may be possible - the question is to find a criterion to decide 
whether it is"^. 

One special case where causality is realised in a particularly simple way is in 
globally hyperbolic spacetimes, where the manifold admits a foliation into a set of 
spacclikc Cauchy surfaces with fibres given by timelike geodesies. It is not hard to 
imagine that the same structure could be preserved using the effective metric G^i, 
to define 'spacelike' or 'timelike', especially if Gj^y is only perturbatively different 
from the actual spacetime metric g^^. But this would be a global question and the 
preservation of global hyperbolicity is not a priori guaranteed. 

The clearest criterion for causality in general involves the concept of stable 
causality discussed, for example, in the monograph of Hawking and EUis"^^. Propo- 
sition 6.4.9 states the required definition and theorem: 

• A spacetime manifold {M,g^,y) is stably causal if the metric ^^u^ has an open 
neighbourhood such that A4 has no closed timclike or null curves with respect to 
any metric belonging to that neighbourhood. 

• Stable causality holds everywhere on M if and only if there is a globally defined 
function / whose gradient D^f is everywhere non-zero and timelike with respect 

to gfj,^. 

'^Bimetric theories of gravity have an extensive literature. See ref.^* for an elegant recent con- 
struction and references therein for earlier work. 

"^The equivalent discussion of causality for superluminal propagation in Minkowski spacetime with 
Casimir plates is given in ref.^^. This also provides a nice example of the distinction between ray 
and phase velocities discussed above. 
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According to this theorem, the absence of causahty violation in the form of 
closed timelike or lightlike curves is assured if we can find a globally defined func- 
tion / whose gradient is timclikc with respect to the effective metric G^i, for light 
propagation. / then acts as a global time coordinate. 

To see how this criterion can be applied to a particular example, one for which 
stable causality is preserved by the new light cone metric, consider the cosmological 
Friedmann-Robertson- Walker spacetime. Since the FRW metric is Weyl flat, the 
modified light cone condition Eq.(3) reads simply 



where C = ^ (2& + 4c) and the energy- momentum tensor is 

= (p + P)nij,nu - PQixu 



(18) 



(19) 



with specifying the time direction in a comoving orthonormal frame, p is the 

energy density and P is the pressure, which in a radiation-dominated era arc related 
by /? — 2>P = 0. The phase velocity is independent of polarisation and is found to 
be super luminaP: 



^ph = ^ = l + ic(P + P) 

At first sight, this looks surprising given that fc^ 
energy condition T^^fe'^fc'^ > 0. However 
along the rays, jf" = Q^^k^, we find 

and 



(21) 



> 0, its sign fixed by the weak 
if instead we consider the momentum 



ray 



The effective metric G = 



G 



/1 + Cp 




\ 



^ is (in the orthonormal frame) 




-(1-CP) 









-(1-CP) 





(22) 



(23) 



(24) 



-{I -CP) J 



and 



In this case, therefore, we find equal and superluminal velocities tiph 
< is manifestly spacelike as required. 

Is stable causality preserved? In this case the answer is yes, since we may still 
use the cosmological time coordinate t as the globally defined function /. We need 



'^In the radiation dominated era, where p{t) = oo— * we have 



Vph 



1+ — C t" 
167r 



(20) 



which, as already observed in ref.^, increases towards the early universe. Although this expression 
is only reliable in the perturbative regime where the correction term is small, it is intriguing that 
QED predicts a rise in the speed of light in the early universe. It is interesting to speculate whether 
this superluminal effect persists for high curvatures near the initial singularity and whether it could 
play a role in resolving the horizon problem in cosmology 
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only check that D^t defines a timehke vector with respect to the effective metric 
G^i,. This is true provided Goo > 0, which is certainly satisfied by Eq.(24). So at 
least in this case, superluminal propagation is compatible with causality. 

2.5 Time Machines? 

Although we have seen that causality is not necessarily violated by superluminal 
propagation, it is important to look for counter-examples where the Drummond- 
Hathrell effect may create a time machine. One imaginitive suggestion was put 
forward by Dolgov and Novikov (DN)^-'' , involving two gravitating sources in relative 
motion. This scenario therefore has some echoes of the Gott cosmic string time 
machine^^; both are reviewed in ref.^^. 

The DN proposal is to consider first a gravitating source with a superluminal 
photon following a trajectory which we may take to be radial. Along this trajectory, 
the metric interval is effectively two-dimensional and DN consider the form 

ds"^ = A'^{r)dt'^ - B^{r)dr^ (25) 

(An explicit realisation is given by radial superluminal signals in the Bondi-Sachs 
spacetime, described in ref.''.) The photon velocity in the {t, r) coordinates is taken 
to be w = 1 + Sv, so the effective light cones lie pcrturbatively close to the geometric 
ones. The trajectory is forward in time with respect to t. 

DN now make a coordinate transformation corresponding to a frame in relative 
motion to the gravitating source, rewriting the metric interval along the trajectory 
as 

ds^ = A^{t',r'){dt'^ -dr'^) (26) 



The transformation is-'^ 



'y{u)(t — ur — uf{r)) 
= 'y{u){r-ut + f{r)) 



(27) 



with 



fir) = y dr (I - l) (28) 
Now, a superluminal signal with velocity 

Bdr 

v = l + Sv=-- (29) 

emitted at (ti,ri) and received at (^2,^2) travels forward in t time (for small, 
positive Sv) with interval 



t2 



-ti= jj dr {I- 5v) J (30) 



■'^This transformation comprises two steps. First, since any 2-dim metric is conformally flat, we 
can bring the metric Eq.(25) into standard form ds^ = (^dP —df^) . Then, a boost with velocity 
u is made on the flat coordinates {t,f) to give the DN coordinates {t',r'). 
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X 



Y 



Figure 3. The Dolgov-Novikov time machine proposal. A superhiminal signal from X, described 
as backwards-in-time in a relevant frame, is sent towards a second gravitating source Y moving 
relative to X and returned symmetrically. 

As DN show, however, this motion is backwards in t' time for sufficiently large u, 
since the equivalent interval is 



The required frame velocity is m > 1 — 6v, i.e. since dv is small, u > -. 

The situation so far is therefore identical in principle to the discussion of su- 
perluminal propagation illustrated in Fig. 1. In DN coordinates the outward su- 
perluminal signal is certainly propagating backwards in time, but a reverse path 
with the same perturbatively superluminal velocity would necessarily go sufficiently 
forwards in time to arrive back within the future light cone of the emitter. 

At this point, however, DN propose to introduce a second gravitating source 
moving relative to the first, as illustrated in Fig. 3. They now claim that a superlu- 
minal photon emitted with velocity v{r) in the region of X will travel backwards in 
time (according to the physically relevant coordinate t') to a receiver in the region 
of Y. A signal is then returned symmetrically to be received at its original position 
in the vicinity of X, arriving, according to DN, in its past. This would then be 
analogous to the situation illustrated in Fig. 2. 

However, as we emphasised in section 2.1, we are not free to realise the sce- 
nario of Fig. 2 in the gravitational case, because the SEP-violating superluminal 
propagation proposed by Drummond and Hathrell is pre-determined, fixed by the 
local curvature. The t' frame may describe back-in-time motion for the outward 
leg, but it does not follow that the return path is similarly back-in-time in the same 
frame. The appropriate special relativistic analogue is the scenario of Fig. 1, not 
Fig. 2. This critique of the DN time machine proposal has already been made by 
Konstantinov^^ and further discussion of the related effect in flat spacetime with 
Casimir plates is given in ref.^^. The relative motion of the two sources, which at 
first sight seems to allow the backwards-in-time coordinate t' to be relevant and to 
be used symmetrically, does not in fact alleviate the problem. 

The true situation seems rather more to resemble Fig. 4. With the gravitating 
sources X and Y sufficiently distant that spacetime is separated into regions where it 
is permissible to neglect one or the other, a signal sent from the vicinity of X towards 
Y and back would follow the paths shown. But it is clear that this is no more than 
stitching together an outward plus inward leg near source X with an inward plus 
outward leg near Y. Since both of these are future-directed motions, in the sense of 
Fig. 1, their combination cannot produce a causality- violating trajectory. If, on the 
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Figure 4. A decomposition of the paths in Fig. 3 for well-separated sources 



other hand, we consider X and Y to be sufficiently close that this picture breaks 
down, we lose our ability to analyse the Drummond-Hathrell effect, since we would 
need the full collision metric for the gravitating sources which is not known for 
physically realisable examples. 

We therefore conclude that the Dolgov-Novikov time machine docs not work. 
The essential idea of trying to realise the causality-violating special relativistic 
scenario of Fig. 2 by using two gravitational sources in relative motion does not in 
the end succeed, precisely because the physical Drummond-Hathrell light cones are 
fixed by the local curvature. Once more it appears that in general relativity with 
SEP-violating interactions, superluminal photon propagation and causality can be 
compatible. 



2.6 The Event Horizon 

We have seen that when quantum effects are taken into account, the physical light 
cones need not coincide with the geometrical null cones. This immediately raises the 
question of black hole event horizons - do the physical horizons for light propagation 
also differ from the geometrical horizons, and are they polarisation dependent? If so, 
this would have profound repercussions for phenomena such as Hawking radiation. 
The answer is best seen using the Newman-Penrose form of the light cone, viz. 

^^-^^^±^^$00 ± '-^{^. + %) = (32) 

If we define the tetrad with as an outward-directed null vector orthogonal to 
the horizon 2-surface, then a fundamental theorem on horizons states that both 
$00 and ^'o are zero precisely at the horizon. The detailed proof, which is given in 
ref.^^, involves following the convergence and shear of the generators of the horizon. 
In physical terms, however, it is easily understood as the requirement that the flow 
of both matter (given by the Ricci term) and gravitational radiation (given by the 
Weyl term) are zero across the horizon. 

It follows that for outward-directed photons with fc^ = ^^^J.^ the quantum cor- 
rections vanish at the horizon and the light cone coincides with the null cone. 
The geometrical event horizon is indeed the true horizon for physical photon 
propagation'*'^'^. Again, no conflict arises between superluminal propagation and 
essential causal properties of spacetime. 
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3 Causality, Characteristics and the 'Speeds of Light' 

So far, our analysis of photon propagation has been based entirely on the leading- 
order, Drummond-Hathrell effective action Eq.(l). However, as we show in section 
4, the full effective action contains terms to all orders in a derivative expansion and 
these must be taken into account to go beyond the low-frequency approximation. 
Photon propagation in QED in curved spacetime is therefore dispersive and we 
must understand how to identiiy the 'speed of light' which is relevant for causality. 



3.1 'SpeeAs of Light' 

An illuminating discussion of wave propagation in a simple dispersive medium is 
given in the classic work by Brillouin^^. This considers propagation of a sharp- 
fronted pulse of waves in a medium with a single absorption band, with refractive 
index n(a;): 

n\u^) = 1 - ^ (33) 

where a, p are constants and cjq is the characteristic frequency of the medium. Five^ 
distinct velocities are identified: the phase velocity v^h. = = ^TJ^' group velocity 
Vgp = -j^, signal velocity Vsig, energy-transfer velocity Vcn and wavefront velocity 
Vvfi, with precise definitions related to the behaviour of contours and saddle points in 
the relevant Fourier integrals in the complex w-plane. Their frequency dependence 
is illustrated in Fig. 5. As the pulse propagates, the first disturbances to arrive 
are very small amplitude waves, 'frontrunners', which define the wavefront velocity 
v^f. These are followed continuously by waves with amplitudes comparable to the 
initial pulse; the arrival of this part of the complete waveform is identified in rcf.^-'^ 
as the signal velocity Vsig- As can be seen from Fig. 5, it essentially coincides with 
the more familiar group velocity for frequencies far from u)q, but gives a much more 
intuitively reasonable sense of the propagation of a signal than the group velocity, 
whose behaviour in the vicinity of an absorption band is relatively eccentric.'' As 
the figure makes clear, the phase velocity itself also does not represent a 'speed of 
light' relevant for considerations of signal propagation or causality. 

The appropriate velocity to define light cones and causality is in fact the wave- 
front velocity Wwf • (Notice that in Fig. 5, ?;wf is a constant, equal to c, independent 
of the frequency or details of the absorption band.) This is determined by the 
boundary between the regions of zero and non-zero disturbance (more generally, a 
discontinuity in the first or higher derivative of the field) as the pulse propagates. 
Mathematically, this definition of wavefront is identified with the characteristics of 
the partial differential equation governing the wave propagation^^. Our problem is 
therefore to determine the velocity associated with the characteristics of the wave 



3In fact, if we take into account the distinction discussed in section 2 between the phase velocity 
Vpii and the my velocity Vray, and include the fundamental speed of light constant c from the 
Lorentz transformations, we arrive at seven distinct definitions of 'speed of light'. 
''Notice that it is the group velocity which is measured in quantum optics experiments which find 
light speeds of essentially zero^^ or many times c A particularly clear description in terms of 
the effective refractive index is given in ref.^^. 
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Figure 5. Sketch of the behaviour of the phase, group and signal velocities with frequency in 
the model described by the refractive index Eq. (33). The energy-transfer velocity (not shown) is 
always less than c and becomes small near loo- The wavefront speed is identically equal to c. 

operator derived from the modified Maxwell equations of motion appropriate to the 
new effective action. 

Notice that a very complete and rigorous discussion of the wave equation in 
curved spacetime has already been given in the monograph by Friedlander^^ , in 
which it is proved (Theorem 3.2.1) that the characteristics are simply the null 
hypersurfaces of the spacetime manifold, in other words that the wavefront always 
propagates with the fundamental speed c. However, this discussion assumes the 
standard form of the (gauge-fixed) Maxwell wave equation (cf. ref.^^, eq. (3.2.1)) and 
does not cover the modified wave equation derived from the action Eq.l, precisely 
because of the extra curvature couplings which lead to the effective metric G^^ and 
superluminal propagation. 

3.2 Characteristics, Wavefronts and the Phase Velocity Wph(oo) 

Instead, the key result which allows a derivation of the wavefront velocity is derived 
by Leontovich^^. In this paper', an elegant proof is presented for a very general set 
of PDEs that the wavefront velocity associated with the characteristics is identical 
to the w — > CX3 limit of the phase velocity, i.e. 

Uwf = hm = lim Vph{uj) (34) 

The proof is rather formal, but is of sufficient generality to apply to our discussion of 
photon propagation using the modified effective action of section 4. We reproduce 
the essential details below. 

*I am very grateful to A. Dolgov, V. Khoze and I. Khriplovich for their help in obtaining and 
interpreting ref.^^. 
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The first step is to recognise that any second order PDE can be written as 
a system of first order PDEs by considering the first derivatives of the field as 
independent variables. Thus, if for simplicity wc consider a general second order 
wave equation for a field u{t, x) in one space dimension, the system of PDEs we 
need to solve is 







(35) 



where .^^ = {«, f , 



at ' Sec J 

Making the 'geometric optics' ansatz 



(j}i = (pi exp i{wt — kx) (36) 

where the frequency-dependent phase velocity is Vph{k) = u}{k)/k, and substituting 
into Eq.(35) we find 

(iuiaij — ikbij + Cij^ipj = (37) 
The condition for a solution. 



det 



aijVpi,{k) - bi 







(38) 



then determines the phase velocity. 

On the other hand, wc need to find the characteristics of Eq.(35), i.e. curves C 
on which Cauchy's theorem breaks down and the evolution is not uniquely deter- 
mined by the initial data on C. The derivatives of the field may be discontinuous 
across the characteristics and these curves are associated with the wavefronts for 
the propagation of a sharp-fronted pulse. The corresponding light rays are the 
'bicharacteristics'. (See, for example, ref.^^ chapters 5.1, 6.1 for further discussion.) 

Wc therefore consider a characteristic curve C in the {t, x) plane separating 
regions where (pi = (ahead of the wavefront) from ^ (behind the wave- 
front). At a fixed point {to,xo) on C, the absolute derivative of (pi along the curve, 
parametrised as x{t), is just 



dt 



dt 



dx 



dx 
dt 



(39) 



where dx/dt = v^,! gives the wavefront velocity. Using this to eliminate from 
the PDE Eq.(35) at {to,xo), we find 



/ dx , \ d(j)j 



4»> 







dt ■ ^ (40) 

Now since C is a wavefront, on one side of which (pi vanishes identically, the second 
two terms above must be zero. The condition for the remaining equation to have a 
solution is simply 



det 



= 



(41) 



which determines the wavefront velocity Vy,f. The proof is now evident. Comparing 
Eqs.(38) and (41), we clearly identify 



Vwf = Vph{k oo) 



(42) 
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The wavefront velocity in a gravitational background is therefore not given a 
priori by c. Taking vacuum polarisation into account, there is no simple non- 
dispcrsivc medium corresponding to the vacuum of classical Maxwell theory in 
which the phase velocity represents a true speed of propagation; for QED in curved 
spacetime, even the vacuum is dispersive. In order to discuss causality, we therefore 
have to extend the original Drummond-HathrcU results for Wph(ti' ^ 0) to the high 
frequency limit Vp\^{uj oo), as already emphasised in their original work. 

A further subtle question arises if we write the standard dispersion relation for 
the refractive index n{ui) in the limit w ^ oo: 



For a conventional dispersive medium, 3n(a;) > 0, which implies that n(oo) < 
n(0), or equivalently Wph(oo) > Wph(O). Evidently this is satisfied by Fig. 5. The 
key question though is whether the usual assumption of positivity of holds 
in the present situation of the QED vacuum in a gravitational field. If so, then 
(as already noted in ref.^) the superluminal Drummond-Hathrell results for t;ph(0) 
would actually be lower bounds on the all-important wavefront velocity Wph(oo). 
However, it is not clear that positivity of 3n(a;) holds in the gravitational context. 
Indeed it has been explicitly criticised by Dolgov and Khriplovich in refs.^'''^^, who 
point out that since gravity is an inhomogcncous medium in which beam focusing as 
well as diverging can happen, a growth in amplitude corresponding to 3n(a;) < 
is possible. The possibility of Wph(oo) < Wph(O), and in particular t)ph(oo) = c, 
therefore does not seem to be convincingly ruled out by the dispersion relation 



4 Dispersion 

After these general considerations, we now return to QED in curved spacetime and 
use the full eff'ective action to study dispersion and investigate the high-frequency 
limit. 

4-1 Effective Action for Photon- Gravity Interactions 

The local effective action in QED to 0{RFF) keeping terms of all orders in deriva- 
tives is derived in ref.®. The result is: 




(43) 



Eq.(43). 




(44) 
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This was found by adapting a background field action valid to third order in gen- 
eralised curvatures due to Barvinsky, Gusev, Zhytnikov and Vilkovisky^^ (see also 
rcf.''") and involves rc-cxprcssing their more general result in manifestly local form 
by an appropriate choice of basis operators. 

In this formula, the G„ (n > 1) are form factor functions of three operators: 



r,2 r)2 n2 

G:^gJ%%^) (45) 
V m'' ) 



where the first entry (£^^i)) acts on the first following term (the curvature), etc. Go 
is similarly defined as a single variable function. These form factors arc found using 
heat kernel methods and are given by 'proper time' integrals of known algebraic 
functions. Their explicit expressions can be found in ref.^. Evidently, Eq.(44) re- 
duces to the Drummond-Hathrell action if we neglect all the higher order derivative 
terms. 



^.2 Dispersion and the Light Cone 

The next step is to derive the equation of motion analogous to Eq.(ll) from this 
generalised effective action and to apply geometric optics to find the corresponding 
light cone. This requires a very careful analysis of the relative orders of magnitudes 
of the various contributions to the equation of motion arising when the factors of 
D'^ in the form factors act on the terms of 0{RF). These subtleties are explained 
in detail in ref.^. The final result for the new effective light cone has the form 

- -^F('^)R^xk^k^ + l-G(^)R^,xpk^k^a''a'' = (46) 

where F and G arc known functions with well- understood asymptotic properties^. 
Clearly, for agreement with Eq.(12), we have F(0) = 26, G(0) = 8c. 

The novel feature of this new light cone condition is that F and G are functions 
of the operator k.D acting on the Ricci and Riemann tensors.-' So although the 
asymptotic behaviour of F and G as functions is known, this information is not 
really useful unless the relevant curvatures are eigenvalues of the operator. On the 
positive side, however, k.D does have a clear geometrical interpretation — it simply 
describes the variation along a null geodesic with tangent vector fc'^. 

The utility of this light cone condition therefore seems to hinge on what we 
know about the variations along null geodesies of the Ricci and Riemann (or Weyl) 
tensors. It may therefore be useful to re-express Eq.(46) in Newman-Penrose form: 

where F = 2F + G. 

Unfortunately, we have been unable to find any results in the relativity literature 
for l.D^Qo and i.D^o which are valid in a general spacetime. In particular, this 

■'Note that because these corrections are already of 0(a), we can freely use the usual Maxwell 
relations k.Dk" = and k.Da" = in these terms; we need only consider the effect of the operator 
k.D acting on R^i, and R^xvXp- 
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is not one of the combinations that are constrained by the Bianchi identities in 
Newman-Penrose form (as displayed for example in ref.^^, chapter 1, Eq.(321)). To 
try to build some intuition, wc have therefore fookcd at particular cases. The most 
interesting is the example of photon propagation in the Bondi-Sachs metric^^'^'^ 
which we recently studied in detail^. 

The Bondi-Sachs metric describes the gravitational radiation from an isolated 
source. The metric is 

ds^ = -Wdu^ - 2e^^dudr + r^hij{dx' - U'du){dx' - U'du) (48) 

where 

hijdx'dx^ = i(e2T+e2'5)d6'2+2sinh(7-^)sin6'd6'd(/>+i(e-2T+e-2'5)sin2 6*^02 (49) 

The metric is valid in the vicinity of future null infinity X+ . The family of hyper- 
surfaces u = const are null, i.e. g'^'^dfj^ud^u — 0. Their normal vector £^ satisfies 

= d^u ^ e^ = o, e^Df^r = o (50) 

The curves with tangent vector £^ are therefore null geodesies; the coordinate r is a 
radial parameter along these rays and is identified as the luminosity distance. The 
six independent functions VF, /3, 7, (5, t/* characterising the metric have expansions 
in i in the asymptotic region near , the coefficients of which describe the various 
features of the gravitational radiation. 

In the low frequency limit, the light cone is given by the simple formula Eq.(4) 
with <l>oo = 0. The velocity shift is quite different for the case of outgoing and 
incoming photons'". For outgoing photons, = i^, and the light cone is 

. ^ 0(1) ,51) 

while for incoming photons, k^ = n^, 

= ±^ (*. + *:) ~ 0(1) (52) 

Now, it is a special feature of the Bondi-Sachs spacetime that the absolute 
derivatives of each of the Weyl scalars ^'o, • • • , 4*4 along the ray direction £^ van- 
ishes, i.e. ^'o,-.-,4'4 are parallel transported along the rays'^"'''^"'. In this case, 
therefore, we have: 

e-D^o = e-D^4 = (53) 

but there is no equivalent simple result for either n ■ D ^'4 or n ■ D ^o- 

Although it is just a special case, Eq.(53) nevertheless leads to a remarkable 
conclusion. The full light cone condition Eq.(47) applied to outgoing photons in 
the Bondi-Sachs spacetime now reduces to 

± £^G(0)(*o + *S) = (54) 

since i-D'^Q = 0. In other words, the low-frequency Drummond-Hathrell prediction 
of a superluminal phase velocity Uph(O) is exact for all firequencies. There is no 
dispersion, and the wavefront velocity t;ph(oo) is indeed superluminal. 
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This is potentially a very important result. Based on the improved effective 
action Eq.(44), we have shown there is at least one example in which the wavefront 
truly propagates with supcrluminal velocity. Quantum effects have indeed shifted 
the light cone into the geometrically spacelike region. 



4-3 Non-Perturbative Effective Action and High- Frequency Propagation 

Unfortunately, there is one final twist to the story which could invalidate the above 
conclusion. If instead of a gravitational field, we consider photon propagation in a 
constant background magnetic field, we find the following birefringent modification 
to the light cones:^^'^*''^ 

e - ^C^Y f durdssN^\,^{u,z)eM^+^'^'nu,.)) ^ (55) 

ZTT V m / J _i Jq 

where z = and Q = ^ — . Exact expressions are known for the functions N 
and P as well as their weak-field expansions. 

Now, in the weak-field, low-frequency regime, Q. ^1 and we can simply expand 
the exponential to leading order in Q^P. The low-frequency limit of the phase 
velocity for the two polarisations is the well-known result^'^''' 

a /eB\2 rl4 8 1 

and is in each case subluminal. At high frequencies, however, Eq.(55) is domi- 
nated by the rapidly-varying phase factor exp(— is'^fJ^P) and the correction to the 
phase velocity tends to zero (from a superluminal value) with a non-analytic 17"^/^ 
behaviour: 

t;ph(a;^oo) = 1 + — j [c||,c^]f]3 (57) 

where C||,cx are known positive constants*. The phase velocity therefore has the 
standard form for a dispersive medium illustrated in Fig. 5. In particular, the 
wavefront velocity v^f = Upii(oo) = c. 

The important lesson for the gravitational case is this. If we had simply used an 
effective action for QED in a background electromagnetic field keeping terms up to 
0{F^) in the field strength and to all orders in derivatives, generalising the Euler- 
Heisenberg Lagrangian, we would have accurately found the leading low-frequency 
dependence of Wph('^) but would have completely missed the non-analytic high- 
frequency behaviour. This arises from terms in the full effective action which are 
non-perturbative in the background field and give rise to the phase factor in Eq.(55). 

If the gravitational case is similar, this would imply that the modified light cone 
can be written heuristically as 

fc2 + « g-»(i+.=oM«,fl)) = (58) 

where both J\f and V can be expanded in powers of curvature, and derivatives of 
curvature, presumably associated with factors of co as in the last section. The fre- 
quency dependent factor Cl would be Q ~ ~ 0{jf^), where 'i?' denotes some 
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generic curvature component and L is the typical curvature scale. If this is true, 
then an expansion of the effective action to 0{RFF), even including higher deriva- 
tives, would not be sufficient to reproduce the full, non-pcrturbativc contribution 
exp(— is'^fi^'P). The Drummond-Hathrell action would correspond to the leading 
order term in the expansion of Eq.(58) in powers of neglecting derivatives, while 
our improved effective action of section 4.1 sums up all orders in derivatives while 
retaining the restriction to leading order in curvature. 

The omission of the non-perturbative contribution would be justified only in 
the limit of small 11, i.e. for <C 1. Neglecting this therefore prevents us from 
accessing the genuinely high frequency limit A ^ needed to find the asymptotic 
limit 'yph(oo) of the phase velocity. Moreover, assuming Eq.(58) is indeed on the 
right lines, it also seems inevitable that for high frequencies (large J7) the rapid 
phase variation in the exponent will drive the entire heat kernel integral to zero, 
ensuring the wavefront velocity v^jf = c. 

4-4 Outlook 

At present, it is not clear how to make further progress. The quantum field theoretic 

calculation required to find such non-pcrturbative contributions to the effective ac- 
tion and confirm an exp(— is'^fi^T-') structure in Eq.(58) appears difficult, although 
some technical progress in this area has been made recently in ref.^* and work in 
progress^''. One of the main difficulties is that since a super luminal effect requires 
some anisotropy in the curvature, it is not sufficient just to consider constant cur- 
vature spacetimes. (Recall that the Ricci scalar term in the effective action Eq.(l) 
docs not contribute to the modified light cone Eq.(2).) A possible approach to 
this problem, which would help to control the plethora of indices associated with 
the curvatures, might be to reformulate the heat kernel calculations directly in the 
Newman-Pcnrosc basis. On the other hand, perhaps a less ambitious goal would be 
to try to determine just the asymptotic form of the non-perturbative contribution 
in the f2 ^ oo limit. 

A final resolution of the dispersion problem for QED in curved spacetime has 
therefore still to be found. If the perturbative expansion of the effective action 
to 0{RFF) is sufficient, then as we have seen there exist at least some examples 
where the wavefront velocity is really superluminal. In this case, all the issues 
concerning causality discussed in section 2 would apply to QED. Perhaps more 
likely, however, is the scenario described above, where the high-frequency dispersion 
is driven by non-perturbative contributions to the effective action such that the 
wavefront velocity remains precisely c. It would then be interesting to see exactly 
how the phase velocity behaves as a function of w (c.f. Fig. 5) and whether the 
explanation advanced in section 3.2 for the non- validity of the standard refractive 
index dispersion relation Eq.(43) is correct. 

Finally, of course, even if it does turn out that Wwf = c for QED itself, the 
discussion of causality in this paper may still be relevant to photon propagation in 
more speculative theories, including bimetric theories of gravity^'*, string-inspired 
non-hnear electrodynamics^", Lorentz and CPT violating effective Lagrangians^" 
and non-commutative gauge theories. 



Shore: submitted to World Scientific on February 7, 2008 20 



Acknowledgments 



This work is supported in part by PPARC grant PP/G/O/2000/00448. 



References 

1. IT. Dnimmond and S. HathrclL Phys. Rev. D22 (1980) 343. 

2. R.D. Daniels and G.M. Shore, Nucl. Phys. B425 (1994) 634. 

3. J. I. Latorre, P. Pascual and R. Tarrach, Nucl. Phys. B437 (1995) 60. 

4. G.M. Shore, Nucl. Phys. B460 (1996) 379. 

5. W. Dittrich and H. Gies, Phys. Lett. B431 (1998) 420-429; Phys. Rev. D58 
(1998) 025004. 

6. R.D. Daniels and G.M. Shore, Phys. Lett. B367 (1996) 75. 

7. G.M. Shore, Nucl. Phys. B605 (2001) 455. 

8. G.M. Shore, Nucl. Phys. B633 (2002) 271. 

9. G.M. Shore, Nucl. Phys. B646 (2002) 281. 

10. K. Scharnhorst, Phys. Lett. B236 (1990) 354. 

11. G. Barton, Phys. Lett. B237 (1990) 559. 

12. S. Liberati, S. Sonego and M. Visser, Annals Phys. 298 (2002) 167. 

13. S.W. Hawking and G.F.R. Ellis, The Large Scale Structure of Spacetime, Cam- 
bridge University Press, 1973. 

14. I.T. Drummond, Phys. Rev. D63 (2001) 043503. 

15. A.D. Dolgov and l.D. Novikov, Phys. Lett. B442 (1998) 82. 

16. J.R. Gott, Astrophys. J. 288 (1985) 422. 

17. G.M. Shore, gr-qc/0210048, to be published in Int. J. Mod. Phys. A 

18. M.Yu. Konstantinov, gr-qc/9810019; Russ. Phys. J. 45 (2002) 23. 

19. S.W. Hawking, 'T/ie Event Horizon', 1972 Les Houches lectures, ed. B. De 
Witt, Gordon and Breach, 1972. 

20. G.W. Gibbons and C.A.R. Herdeiro, Phys. Rev. D63 (2001) 064006. 

21. L. Brillouin, Wave Propagation and Group Velocity, Academic Press (London) 
1960. 

22. L.V. Hau, Scientific American 7 (2001) 66. 

23. L.J. Wang, A. Kuzmich and A. Dogoriu, Nature 406 (2000) 277. 

24. R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol II, Inter- 
science, New York, 1962. 

25. F.G. Priedlander, The Wave Equation on a Curved Spacetime, Cambridge Uni- 
versity Press, 1975. 

26. M.A. Leontovich, in L.I. Mandelshtam, Lectures in Optics, Relativity and 
Quantum Mechanics, Nauka, Moscow 1972 {in Russian). 

27. A.D. Dolgov and LB. Khriplovich, Sov. Phys. JETP 58(4) (1983) 671. 

28. LB. Khriplovich, Phys. Lett. B346 (1995) 251. 

29. A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky and V.V. Zhytnikov, Print-93- 
0274 (Manitoba), 1993. 

30. A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky and V.V. Zhytnikov, J. Math. 
Phys. 35 (1994) 3525; J. Math. Phys. 35 (1994) 3543; Nucl. Phys. B439 
(1995) 561. 



Shore: submitted to World Scientific on February 7, 2008 21 



31. S. Chandresekhar, The Mathematical Theory of Black Holes, Clarendon, Ox- 
ford, 1985. 

32. H. Bondi, M.G.J, van der Burg and A.W.K. Metzner, Proc. Roy. Soc. A269 
(1962) 21 

33. R.K. Sachs, Proc. Roy. Soc. A270 (1962) 103. 

34. R.A. d'Inverno, Introducing Einstein's Relativity, Clarendon, Oxford, 1992. 

35. W. Tsai and T. Erber, Phys. Rev. DIO (1974) 492. 

36. W. Tsai and T. Erber, Phys. Rev. D12 (1975) 1132. 

37. S. Adler, Ann. Phys. (N.Y.) 67 (1971) 599. 

38. A.O. Barvinsky and V.F. Mukhanov, Phys. Rev. D66 (2002) 065007. 

39. Yu.V. Gusev, private communication. 

40. V.A. Kostelecky, Phys. Rev. D58 (1998) 116002. 



Shore: submitted to World Scientific on February 7, 2008 22 



